Spin-Peierls instabilities of antiferromagnetic rings in a magnetic field 



Valeria Lantc,^'^^* loannis Rousochatzakis,^ Karlo Penc,'^ Oliver Waldmann,'' and Frederic Mila^ 

^ Dipartimento di Fisica e Matematica, Universita deU'Insuhna, Via Valleggio 11, 1-22100 Como, Italy 
^Institut de theorie des phenomenes physiques, Ecole Polytechnique Federale de Lausanne, CH-1015 Lausanne, Switzerland 
'^Research Institute for Solid State Physics and Optics, H-1525 Budapest, P.O. Box 49, Hungary 
Physikalisches Institut, Universitdt Freiburg, Hermann- Herder-Strafie 3, 79104 Freiburg, Germany 

(Dated: December 18, 2008) 

Motivated by the intriguing properties of magnetic molecular wheels at field induced level cross- 
ings, we investigate the spin-Peierls instability of antiferromagnetic rings in a field by exact diagonal- 
izations of a microscopic spin model coupled to the lattice via a distortion dependent Dzyaloshinskii- 
Moriya interaction. We show that, beyond the unconditional instability at level crossings for in- 
finitesimal magnetoelastic coupling, the model is characterized by a stronger tendency to distort 
at higher level crossings, and by a dramatic angular dependence with very sharp torque anomalies 
when the field is almost in the plane of the ring. These predictions are shown to compare remarkably 
well with available torque and Nuclear Magnetic Resonance data on CsFes. 

PACS numbers: 75.10.Jm, 71.27.+a, 74.20.Mn 



Intermediate between single spins and bulk magnets, 
molecular magnets have attracted a lot of attention since 
they offer a suitable platform for probing the predictions 
of quantum mechanics, e.g. for the tunneling probability 
between almost "classical" states[l]. At the field-induced 
level crossings (LCs), a small gap is in most cases opened 
by small anisotropics, allowing the system upon sweeping 
a magnetic field to remain in the same state or to tunnel 
following the adiabatic ground state [1, 2, 3, 4]. The sub- 
class of ring-like molecular magnets is exceptional in that 
respect. The structure of these "magnetic wheels" is in 
general so symmetric that, although present, anisotropic 
interactions arc predicted to leave intact the degener- 
acy at LCs. However, this prediction is in contradic- 
tion with several experiments. The direct evidence of 
level repulsion in Fe6:Li clusters lead Affronte et aL[5] 
to postulate a distortion at low temperatures that would 
allow extra terms such as Dzyaloshinskii-Moriya (DM) 
interactions [6]. Cinti et al.[7] have introduced a model 
with rigid dimerization and fixed DM anisotropy to ac- 
count for the tunnel splittings in Fe6:Na clusters, while 
Miyashita et al.[8] have discussed the possibility of DM 
anisotropy induced by thermal fluctuations. 

More recently, for the octanuclear CsFeS cluster, a pe- 
culiar torque signal was reported[9], that appears quite 
abruptly around the lowest LC field. This behaviour is 
clearly inconsistent with the previous mechanisms, and 
the scenario of a field-induced magneto-elastic instability 
was suggested[9]. In the same field regime, Schnelzer et 
aL[10] found a large broadening of the proton Nuclear 
Magnetic Resonance (NMR) spectrum which signals the 
presence of large staggered transverse polarizations [10], 
an effect which points to DM anisotropy[ll]. On the 
theoretical side, the experiments were first interpreted in 
the context of a phenomenological 2-level approach[9, 12] 
which results, at the level of a single molecule, in a struc- 
tural instability around the LC due to an anisotropy in- 



duced off-diagonal coupling between the two lowest spin 
levels. This model successfully reproduced the experi- 
mental findings qualitatively, but key questions such as 
which type of magnetic anisotropy triggers the instabil- 
ity are beyond its reach. More recently, Soncini et al.[13\ 
discussed a model which also relies on phenomenological 
off-diagonal coupling parameters, but is based on elastic 
interactions between molecules and attributes the insta- 
bility at the LC to the on-site magnetic anisotropy. The 
semi-empirical nature of these models calls for a micro- 
scopic approach which should also answer two fundamen- 
tal questions: i) does a spontaneous structural instabil- 
ity occur also in a microscopic model, and ii) what is 
the role of (staggered) DM interactions, which should be 
expected to provide stronger off-diagonal couplings than 
on-site anisotropy as they couple directly the two spin 
levels. 




(a) dimerization 



(b) bucl^ling 



FIG. 1: (color on-line) Pictorial representation of the dimer- 
ized (a) and the buckled octagon (b). The d vectors are de- 
termined by the symmetry generators. In (a) these are the C4 
axis, the vertical refiection plane (t„ and the horizontal plane 
ah (here the generated group is = D4 x i), while in (b) 
they are the rotorefiection Sg = fhCg and the C2 axis. 
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In this Letter, we investigate the possibihty of a struc- 
tural instabiUty in the context of a fully microscopic 
model where a general A^-site molecular wheel is de- 
scribed by a spin-s Hamiltonian coupled to the lattice. 
This has allowed us to prove that molecular wheels are 
indeed unconditionally unstable at LCs. For the case 
of CsFcs, the general trends around the first LC show 
a remarkable agreement with all available experimental 
data[9, 10] for fields almost perpendicular to the plane of 
the ring, lending strong support to the present instabil- 
ity scenario. In addition, the model allows us to compare 
what happens at different LCs and field orientations and 
to make a number of predictions which go beyond the pre- 
vious approaches. In particular, we find that the model is 
characterized by a much stronger tendency to distort at 
higher LCs. It also shows a dramatic angular dependence 
whereby the torque shows a step-wise linear behavior or 
sharp anomalies for fields almost perpendicular or paral- 
lel to the ring plane respectively. The latter is found to 
be at least twice larger than what expected from the low- 
est order perturbation theory, showing that sub-leading 
corrections play a substantial role. 

The Hamiltonian we consider is adapted from the mod- 
els used to describe the spin-Peierls transition in spin- 1/2 
chains, with two important differences: (i) We add a sin- 
gle ion anisotropy since we are interested in spins larger 
than 1 /2 (this term is also crucial to reproduce the char- 
acteristic background torque signal of Ref. [9]), and (ii) 
we include a distortion-dependent DM interaction. The 
Hamiltonian reads: 

n = ^j,s, •s,+i-i?^4-B-s 

■i i 

+ yH-^^ +11^, -s, xs,+i (1) 

i i 

with implicit periodic boundary conditions. The nearest- 
neighbor exchange is antiferromagnetic {Ji > 0) and the 
single-ion anisotropy easy-axis {D > 0). Si denotes a 
deformation parameter of the bond (in the simplest case 
it is the variation of the length of the bond) between 
the sites i and i + 1, and di is the corresponding DM 
vector, assumed to depend on Si. In what follows, we 
define a fixed xyz reference frame with the molecule on 
the xy-plane (cf. Fig. 1) and the field B in the xz-plane 
subtending an angle 6 from the z-axis. 

In the spin-Peierls transition of spin-1/2 chains, the 
structure factor diverges at the zone boundary, and an in- 
stability can be triggered by an infinitesimal dimcrization 
of either the exchange integral[14] or the DM interaction 
if it is allowed by symmetry [15, 16]. Since the exchange 
is typically much larger, its effect usually dominates. In 
the present case of a finite system, the situation is quite 
different. In the absence of distortion, the spectrum is 
gapped except at the level crossings, and an infinitesimal 
distortion can lead to a structural instability if and only 



if three conditions are met: (i) One sits at a level cross- 
ing, (ii) the induced perturbation couples the two levels, 
and (iii) the coupling is of first order in the distortion to 
overcome the cost in elastic energy. Now, the invariance 
of the Hamiltonian under cyclic permutations Cn of the 
spin indices implies that the eigenstates can be labeled by 
the momentum k = 2Tm/N (71 = 0, . . . , iV — 1) with the 
ground state alternating from /c = to fc = tt at subse- 
quent LC fields[17] . By symmetry, the last two conditions 
can only be met by a dimerized DM interaction along z 
(sec Fig. 1(a)). Indeed, a modulation of the exchange 
integrals (as studied e.g. in Ref. 18) leads to an SU(2) 
invariant perturbation and cannot lift the degeneracy be- 
tween states with different total spin, while the in-plane 
DM interactions allowed by buckling (see Fig. 1(b)) leads 
to a perturbation Vb which only couples the two levels to 
order 7V/2, and only when N/2 is even[19]. In addition, 
as for the exchange in the spin-Peierls mechanism, the 
dimerization of the DM interaction is expected to be lin- 
ear in 5 in the limit of small S so that condition (iii) is 
met. This leads to the minimal set of conditions: Ji = J, 
Si ~ (—1)^(5, and dj = (— l)*c?((5) e^, with d{S) oc S when 
(5^0. All other effects are irrelevant in the limit of small 
Si and need not be considered further given the purpose 
of the present paper. Our model is further specified by 
the following conventions. All energies arc in units of 
,7=1, while the scale of S is implicitly set by the value 
of K. Besides, we impose a cut-off on the DM interaction 
by choosing d = dotanhJ with do = 0.05 J. Physically, 
such a cut-off must be imposed since DM interactions 
are typically at most a few percent of the exchange. Its 
specific form has been chosen for numerical convenience 
but docs not affect the results qualitatively. 

The model is solved in two steps: (i) Determine the 
ground state energy E{S) for a given value of the dis- 
tortion S. This has been achieved with exact diagonal- 
izations based on a Lanczos algorithm; (ii) Find the dis- 
tortion ^0 which minimizes the energy. This has been 
achieved by a bisection method. We have systematically 
investigated various N, s, B, 9 and K. The main qual- 
itative features are similar for all ring sizes and spins s. 
Here, in order to make contact with the reported experi- 
ments on CsFes, we restrict ourselves to the spin s = 5/2 
octagon case with D = 0.027 (this corresponds to the 
experimental[9] estimates D ~ 0.56 K and J ~ 20.6 K 
for CsFcg) and take 9 = 93.6°, —3.3° as considered in [9]. 
As to the value of K, this is adjusted so that we get an 
approximate agreement with the width of the reported 
torque anomalies at large angles. 

First of all, let us note that the directions of the lo- 
cal magnetizations mi = (s^) and the torque r are fixed 
by symmetry. Indeed, the plane of d and B contain- 
ing the center of the ring (cf. Fig. 1(a)) is a mirror 
plane[20] which, together with the remaining C4 axis, 
gives {mix,miy,mi^) = {m^+i.^, -rm+i^y.mi+i^^), i.e. 
the moments are uniform in the xz-plane and staggered 



s=5/2 octagon at lowest LC field, 6=93.6° 




b 1 I 1 I I I 1 Id 

0.4 0.45 0.5 0.55 0.6 

B 



FIG. 2: (color on-line) Instability at the lowest LC field of 
CsFe8 for 6 = 93.6° and K = 0.03. (a) Optimal values of 
the lowest gap A and the torque signal r with and without 
the instability, and the DM amplitude d = dotanh^o- (b) 
Local moments of any two neighboring sites (here 1 and 2, cf. 
Fig. 1(a)) in the rotated frame x'yz' where e^/ is along B. By 
symmetry, mi^' and rriiz' are uniform while rrny is staggered. 
As expected, nii^i ~ 1/8 after the LC field which corresponds 
to a total polarization of S-i ~ 1 shared by the 8 Fe'^^ ions. 

along the y-axis. Accordingly, r points along the y-axis. 

Now, let us discuss the results we have obtained for the 
first LC at 6* = 93.6°. As can be seen in Fig. 2(a), a finite 
distortion (and thus DM interaction) appears around the 
LC field. It opens a gap A and the jump in the torque 
T is replaced by an almost linear behavior. The local 
magnetizations for (any) two neighboring spins are de- 
picted in Fig. 2(b). Apart from a linear- like uniform 
response in the xz-plane, we also find a large staggered 
magnetization along the y-axis with \miy \ ~ 0.56 at the 
center of the LC. This should be contrasted with the 
phenomenological two-level model of [10] which gives a 
staggered response in the xz-plane (perpendicular to B) 
which is ^ 1.71 at the center of the LC. This difference 
stems from the phase of the off-diagonal coupling, which 
is purely imaginary in our model and was assumed to be 
real in Ref. [10]. A quantitative fit of the NMR data of 
[10] based on our results for the spin polarizations should 
take into account the anisotropic character of the dipolar 
hypcrfinc field, which is clearly beyond the scope of the 
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field interval which includes the lowest three LC points, 
present work. 

We also note that the transition at the first critical 
field is weakly first order in our model, which induces for 
example a small jump in the torque and in vriiy. This 
is not an artifact of the minimization method. Indeed, 
we have checked that the fourth-order derivative of the 
ground state energy functional E(So) is negative in this 
parameter range. This is not generic however, and the 
nature of the phase transitions as a function of B, and 
K will be investigated in detail in a future article. 

We now turn to the remaining LCs at 6* = 93.6°. Fig- 
ure 3 shows the DM amplitude d as well as the torque r 
and the lowest energy gap A, with and without the in- 
stability up to the third LC field. Quite remarkably, the 
tendency to distort becomes much larger at higher LCs, 
actually so large that after the first critical field of the 
second LC the distortion never disappears. This unex- 
pected effect can be traced back to the fact that the local 
spin polarizations in each x s^- term generally grow at 
large magnetizations (or fields), and by also noting that 
subsequent LC points are closer to each other at high 
fields due to the quasi-continuum character of the high 
energy spectrum[21]. This tendency seems to be at vari- 
ance with the assumption of Ref. [13] that the vibronic 
coupling decreases with the field. However, interestingly 
enough, this effect agrees with the general trend of the 
torque data at higher LCs, where the almost linear be- 
havior is smoothed out. 

Finally, let us discuss the angular dependence of the 
torque (cf. Fig. 4). It turns out to be rather dramatic, 
in fact much more dramatic than anticipated by the phe- 
nomenological approach. For the set of parameters used 
at large angles, the instability for 6 = —3.3° at the first 
LC field could barely be detected, and for the purpose of 
illustration, we have chosen a slightly smaller stiffness K 
to make the degeneracy lifting visible. This tendency is 
again in agreement with the experimental data of Ref. [9] , 
where additional sharp torque anomalies can be seen at 
the second and third LCs, but much weaker at the first 
one. When comparing the shape of the anomalies with 
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FIG. 4: (color on-line) Instability for the lowest three LCs of 
CsFes for 9 = 93.6° and K = 0.02. (a) Torque r and lowest 
gap A with and without the instability, (b) DM amplitude d. 



experiments, one should bear in mind that, for such nar- 
row instabihties, inhomogeneities are likely to lead to a 
distribution of critical fields broader than the instability 
itself, leading to a smoothening and broadening of the 
anomalies. A surprising aspect of our results is that while 
the shape of the distortion agrees with the dome- like con- 
tribution predicted by the lowest order perturbation the- 
ory, their magnitude is at least twice as large [22]. This 
finding can be attributed to the sub-leading perturbative 
correction which scales as dcos^?[ll]. Finally, the dome- 
shape anomalies are expected to appear only when the 
step-height of the torque without the instability is suffi- 
ciently small. This explains why the sharp-anomalies are 
experimentally visible only at small angles 9, except for 
the lowest LC field. 

In conclusion, we have investigated the equivalent of 
the spin-Peierls instability in antiferromagnetic rings in 
the context of a microscopic spin model coupled to the 
lattice via a distortion dependent DM interaction. Apart 
from demonstrating that magnetic rings are unstable at 
all LCs for an infinitesimal distortion, our results ex- 
plain the general trends of all currently available data 
for the lowest LC of CsFes at large and small angles, 
lending strong support to the present DM magnetoelas- 
tic model. Our model has also allowed us to make a 
number of specific predictions which go beyond the previ- 
ous phenomenological approaches. For instance, we find 
that the instabilities are much stronger at high magnetic 
fields and may even persist from the second LC onward 
for fields almost perpendicular to the ring plane. We 



also find sharply peaked torque anomalies when the field 
is almost parallel to the molecule in very homogeneous 
samples. It is our hope that the present paper will trigger 
further experimental investigation of these instabilities. 
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Since C'^^^Vb = TVi, for N/2 even or odd, Vt contains 
representations with n — 1,3, . . . , N — 1 for N/2 even, 
and n = 0, 2, . . . , iV - 2 for N/2 odd. Thus a fc = tt 
momentum transfer necessary for the level mixing will 
first appear in A''/2-th order in d and only for A''/2 even. 
Here the mirror operation stands for a reflection through 
the xz-p\a.ne followed by the time reversal operation T. 
For CsFeg, we have checked numerically that the off- 
diagonal coupling Vj between the lowest levels at the 
j-th LC fleld grow with j. For instance, at 6' = —3.3°, 
v/d ~ (0.1565,0.2858,0.4291) at the lowest three LC 
flelds Be ~ (0.734,1.145,1.6373), while at 6* = 93.6°, 
v/d ~ (0.9251, 1.6327, 2.3154) at the corresponding fields 
Be ~ (0.485,1.052,1.598). For a spin-s tetramer, Vj = 



((4s + 1)2 -f)dsme. 
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[22] For instance, for the second and the third LC field at = values (0.0117, 0.0159) (cf. Fig. 4). 

—3.3°, lowest order perturbation theory predicts dmax — 
(0.0045, 0.0067) to be compared with the corresponding 



